Cronin effect and high-/?^ suppression from the 
Color Glass Condensate 1 



Edmond Iancu 

Service de Physique Theorique, CEA/DSM/SPhT, 91191 Gif-sur-Yvette Cedex, France 

Abstract. I give a pedagogical survey of the nuclear collective effects associated with gluon 
saturation and their impact on particle production in high-energy proton (or deuteron)-nucleus 
collisions at RHIC. At central rapidity, the theory predicts a Cronin peak due to independent multiple 
scattering off the valence quarks in the nucleus. At forward rapidities, the peak flattens out and 
disappears very fast, because of the correlations induced through quantum evolution in the nuclear 
gluon distribution at small x. Also, the ratio R p a between the particle yield in proton-nucleus and 
proton-proton collisions is rapidly suppressed when increasing the rapidity, because of saturation 
effects which slow down the evolution of the nucleus compared to that of the proton. This behaviour 
could be responsible for a similar trend observed in the deuteron-nucleus collisions at RHIC. 



I. INTRODUCTION 

When accelerated up to RHIC energies (~ 100 GeV/nucleon), a heavy ion such as a 
gold nucleus is expected to 'evolve' into a high-density form of hadronic matter — the 
"Color Glass Condensate" (CGC) Jlil2,|3] — whose properties are qualitatively different 
from those of 'ordinary' partonic systems at lower energies. This is the matter made of 
the small-*, 'saturated', gluons, and is characterized by an intrinsic momentum scale, 
the saturation momentum Q s [4], which rises rapidly with the energy and the atomic 
number A (Q^(x,A) ~ x~^A s with X ~ 0.3 and 8 ~ 1 /3) and acts as an infrared cutoff 
in momentum integrals involving the gluonic spectrum (see yD for recent reviews). As 
usual, the variable x denotes the longitudinal momentum fraction of the 'interesting' 
gluons — those which participate in the scattering — , and decreases rapidly when 
increasing the energy of the collision. Thus, for sufficiently small x and/or large A, 
Q s (x, A) is a hard scale (Q^ ^> Aq CD ), which allows us to rely on perturbative techniques 
to get insight into the properties of the CGC and of the high energy scattering in QCD. 

Of course, the CGC remains 'virtual' (it is a part of the nuclear wavefunction) be- 
fore any scattering takes place, but one expects it to be 'liberated' in a high-energy 
nucleus-nucleus (or proton-nucleus) collision, and thus to determine the properties of 
the partonic system which is created immediately after the collision 0>|^]. This system 
must be strongly interacting (because of its high energy density and the many interven- 
ing degrees of freedom), and also highly off-equilibrium (in particular, because of the 
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strong anisotropy in the initial distribution of energy and momentum), so it is expected 
to undergo a violent evolution, characterized in the early stages by a rapid longitudi- 
nal expansion and a redistribution of the energy, momentum and other quantum number 
among the various modes, possibly leading to a thermalized state (the "Quark-Gluon 
Plasma") at intermediate stages q], which then further expands and cools down, until 
it hadronizes into a multitude of hadrons (a few thousands) that are eventually captured 
by the detectors at RHIC 10,0]. It appears therefore as a challenge for the theorists to 
imagine observables which could survive (almost) unchanged to such a violent evolu- 
tion, and thus carry out information about the initial conditions (in particular, about the 
color glass). It is furthermore a challenge for the experimentalists to extract and measure 
such observables with a significant accuracy. 

One of the observables which look most promising for a study of the initial conditions 
(since less sensitive to 'final state interactions') is the total multiplicity of produced 
particles 
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where dA^/d 2 p^d?] denotes the spectrum of the produced hadrons (of a given species, 
or summed over several species), p± is the hadron transverse momentum, and r\ is its 
(pseudo)rapidity : 
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with p = Jp\+p 2 z and 6 the angle between the direction of the produced hadron and 
the collision axis. It is an experimental fact that the integral in Eq. © is dominated by 
small momenta. In fact, in lowest order perturbation theory, the spectrum diverges at 
small p^ like l/p\, and the corresponding integral is ill defined. But the non-linear 
effects encoded in the CGC — which, from the point of view of perturbation theory, 
correspond to all-orders resummations of high parton density effects — provide an 
infrared cutoff in the form of the saturation momentum Q s , which at RHIC is estimated 
as Q s ~ 1 GeV. Note that this is marginally a hard scale, which justifies a posteriori the 
use of perturbation theory in the construction of the theory for the CGC O]. 

One sees that the very calculability of the total multiplicity at RHIC within pertur- 
bative QCD is an indirect evidence for the saturation physics leading to the CGC. Note 
that, implicit in these considerations, is the assumption that the spectrum of the pro- 
duced hadrons (as measured in the experiment) can be identified with the spectrum of 
the partons (mostly gluons) liberated in the early stages of the collisions (as computed 
in the CGC). This 'duality' hypothesis is far from being rigorous, but it has some theo- 
retical motivation, and, more importantly, it seems to be roughly consistent with the data 
in various experimental settings. With this assumption, the CGC prediction for the total 
multiplicity dN/drj at RHIC [8] turns out to be quite close to the actual value measured 
in the experiments. In fact, simple calculations within the CGC framework satisfactorily 
describe the global features of the total multiplicity at RHIC, like its dependencies upon 
the energy y/s and the rapidity rj, and also upon the centrality of the collision (or the 
"number of participants") [9]. 



Another interesting prediction of the CGC theory that will be discussed at length in 
what follows refers to the interplay between the energy dependence and the nuclear 
(A) dependence of the gluon spectrum (by which I mean either the unintegrated gluon 
distribution in the nuclear wavefunction, or the spectrum of the gluons produced in a 
heavy ion collision, and which is further assimilated with the spectrum dN/d 2 p±dr} 
of the produced hadrons). Because of coherence effects associated with the non-linear 
gluon dynamics, the gluon distribution of a large nucleus is not simply the incoherent 
sum of the corresponding distributions produced by A separated nucleons. This is true 
not only in the saturation region at low momenta p± ^ Q s (where, as we shall see, the 
gluon distribution scales roughly like A 2 / 3 rather than like A), but also for momenta well 
above Q s , where however the trend of the A-dependence depends upon y_ and get s re- 
versed when increasing the energy : Whereas at low energies, the theory iflL llCA fllL Il2 1 
predicts an enhancement of the A-dependence due to 'higher-twist' rescattering effects 
(which fall off rapidly with increasing p±; see Sect. Ill below), on the other hand, at 
sufficiently large energies, the quantum evolution of the color glass [2] considerably 
modifies the dominant, leading-twist, contribution (this acquires an 'anomalous dimen- 
sion' fld^l [HI H), and thus provides a power-like suppression of the A-dependence, 
which persists up to relatively high pj_ [13] (see also Sect. IV). This suppression is 
further amplified by running coupling effects 111 611 . 

This change of behaviour has dramatic consequences for the spectrum of the produced 
gluons (or hadrons): At low energies, one expects an enhancement at intermediate p± 
("Cronin peak") in the properly normalized production cross-section lfl7L [l8ll . while 
at higher energies, one should rather find a depletion within a wide range of momenta 
("high-p^ suppression") 

When the latter feature has been conceptually first realized lfl9ll . the "high-p^ sup- 
pression" was freshly discovered at RHIC [6] — the high-p^ yield in gold-gold colli- 
sions was found to be almost an order of magnitude lower than expected from jet pro- 
duction arising from incoherent parton-parton scattering [23] — , so it was tempting to 
propose this "quantum saturation" scenario as a possible explanation of the data 11911 . To 
further substantiate this scenario, it is worth mentioning that the onset of the anomalous 
dimension is closely related to another important prediction of the CGC, namely a new 
scaling law for the gluon distribution H13L fl4L fl5L l24fl . which seems to be confirmed by 
the recent discovery of "geometric scaling" in the HERA data at small x < 10~ 2 [25]. 
But in the context of Au+Au collisions at RHIC, the argument in Ref. id suffers from 
a serious drawback: A central-rapidity (rj = 0) jet with p± = 5 — 10 GeV is produced by 
partons with x ~ 10 _1 , which is a too large value of x for the quantum evolution of the 
CGC to play any role! Besides, the observed high-p^ suppression can be also under- 
stood as arising from a completely different physical mechanism: energy loss through 
final-state interactions (or 'jet quenching') Il26ll . 

Nowadays, we know that this last explanation is in fact the correct one. The decisive 
evidence in that sense came from another set of experiments at RHIC in which one of the 
two heavy nuclei beams has been substituted by a deuteron beam, thus eliminating the 
final-state interactions. Fig. [T] exhibits the experimental results obtained by PHENIX 
and STAR for the nuclear modification factor in deuteron-gold collisions (i?dAu) com- 
pared to gold-gold (i?AuAu)- (BRAHMS and PHOBOS have reported similar results; see 
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FIGURE 1. Central rapidity (tj = 0) PHENIX n° and STAR h ± data comparing R dAu to i? AuAu 
12711 . These results, together with similar data from BRAHMS and PHOBOS |27], prove that jet 
quenching in Au+Au collisions at central rapidity is necessarily a final state effect. 



Refs. in iEtIm Here, for generic nuclei A and B which scatter with impact parameter b±, 
the ratio Rab is defined as: 

„ . , v _ 1 d^VA+B/dT?d 2 p ± d 2 ^ 
Rp&{n,PL,bl_) = -fTj 7777 ~TTj 777772 ^27^' 

{Nc n{b)) dN p+p /drid z p ± d 2 b ± 

where (N co u(b±)} denotes the average number of collisions for a centrality bin at the 
given impact parameter 2 . This definition is such that Rab would be equal to one if 
A+B were an incoherent superposition of nucleon-nucleon collisions; converserly, its 
deviation from one is a measure of collective nuclear effects. The data in Fig. [T]clearly 
show enhancement (rather than suppression) in the R&au ra tio at moderately large p±. 
As already mentioned, and will be detailed in Sect. Ill, this Cronin peak can be attributed 
to the multiple interactions suffered by the colliding parton from the deuteron on its way 
through the nucleus (see Fig. |2j). The parton scatters off gluons with a relatively large 
momentum fraction x ~ 0. 1, but whose density is still quite large (hence the importance 
of multiple scattering), because the gold nucleus contains many color sources already 
at low energy: the 3A ~ 600 valence quarks. The results in Fig. [T] demonstrate that the 
strong suppression (by a factor 4 to 5) observed in i?AuAu cannot be attributed to the 
quantum evolution of the gluon distribution in the nuclear wavefunction (otherwise, a 
similar effect would be seen at rj = also in /?dAu)> and thus indirectly confirm its origin 
as a final-state effect. 

More recently, the BRAHMS experiment has presented the first results for 7?dAu 
at 'forward rapidities' (rj = 2 — 3 in the deuteron fragmentation region; see Fig. |3l) 



2 For 'minimum bias' events which are averaged over all values of b±, (N co \\) scales like ABj (Ra+Rb) 2 - 
E.g.: for A+A, (N co n) - A 4 / 3 , while for d+A, (jV co ii) - A 1 / 3 . 
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FIGURE 2. Multiple scattering in a proton-nucleus collision. In the experiments at RHIC, the 
proton is actually replaced by a deuteron. 
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FIGURE 3. Kinematics for gluon production via gluon fusion in proton(j?)-nucleus(A) colli- 
sions (d+Au collisions are similar). With light-cone notations [k^ = (k + ,k~ ,kj_), with k = 
(Ek ± k z )/^/2], the nucleons involved in the collision have purely longitudinal momenta: 
= (0,Pi,0x) and = (.P2,0,0jJ, with P\ = Pi = y/s/2 in the center-of-mass frame, while 
the emerging gluon has = (^ e?J >'^ e _7, !P±J- Energy-momentum conservation implies 

x\ = ky jP\ = (/?^/\A) e 71 an d x 2 = ^2 /Pi = (p±/ \/s)e. r l , where k\ t 2 refer to the two gluons 
which fuse with each other. Thus, larger positive values for r/ correspond to smaller values for 
the longitudinal fraction x\ of the gluon from the nucleus. 

lE^l . which have been soon after confirmed by the other collaborations [7]. The most 
remarkable feature about these new data is that they show a strong dependence upon 
T], and also a general trend with both r\ and the centrality which at a first sight looks 
counterintuitive and surprising: As argued before, the Cronin enhancement is attributed 
to multiple scatterings between a parton from the deuteron and the gluons in the nucleus. 
When increasing r\ , we are probing gluons with smaller values of x (see kinematics in 
the capture to Fig.|3]), so the gluon distribution must increase, and this should favour the 
rescattering effects and thus enhance the Cronin peak. Similarly, more central collisions 
are probing regions with higher density, which should again enhance the Cronin peak. 
But the experimental results show precisely the opposite trends (cf. Fig. |5J ! 

But as counterintuitive as it might seem, this behaviour has been in fact anticipated, 
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FIGURE 4. BRAHMS results for d+Au collisions at 200 GeV/nucleon jH : The ratio R dAu 
for charged hadrons as a function of p± for central and forward pseudorapidities. 
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FIGURE 5. d+Au collisions at BRAHMS fcSD : The ratio 7? C p of yields from collisions of a 
given centrality class (0-20% or 30-50%) to yields from more peripheral collisions (60-80%), 
scaled by the mean number of binary collisions in each sample, as a function of p± for central 
and forward pseudorapidities. 



on the basis of the CGC ideas |Jj& 2fl, 22] : As already mentioned, one effect of the 
quantum evolution towards smaller values of x is to suppress the dependence of the 
high-p^ spectrum upon the atomic number A, or, more generally, upon the number of 
colliding particles in a given centrality bin. What came nevertheless as a surprise is the 
fact that this suppression sets in so fast ! As manifest on Fig. 13 the centrality dependence 
gets reversed already for r\ = 2 (which corresponds to x ~ 10 -3 for p± = 4- GeV), while 
for r\ = 3.2 (x ~ 5 x 10~ 4 ), the data in Fig. HJshow suppression (RdAu < 1) at all the 
measured values of p±. 

In fact, already before the advent of the BRAHMS data, the rapid evolution of the ratio 
Rab with increasing r\ has been observed theoretical ly | 22ll. in a numerical calculation 
based on the Balitsky-Kovchegov (BK) equation lE9ll30ll . In Ref. lE3l . the Cronin peak 
was seen to flatten out and completely disappear after only one unit of rapidity. But 
the mechanism responsible for such a rapid evolution has been fully elucidated only 
subsequently, through an exhaustive analytic study 113 ill , which has also clarified other 
important issues like the respective roles of the proton and the nucleus in the evolution 
with 7], and the effects of the running of the coupling. The main results of this study will 
be briefly described in the remaining part of this review. 



II. THE tfpA-RATIO IN THE CGC 



For simplicity, in what follows I shall discuss proton(p)+nucleus(A) scattering, and 
assume that similar results hold also for the d+Au collisions of interest at RHIC. The 
physically interesting quantity is then the ratio (compare to Eq. ©) 



dAf p A 

R p a(1,P±) = dNpp , (4) 

dl]d 2 p ± d 2 b ± 

between the yield of produced gluons per unit phase-space (i.e., per unit of pseudora- 
pidity, transverse momentum and impact parameter) in pA and pp collisions normalized 
by A 1//3 (since, at a fixed impact parameter, the number of nucleons available for scat- 
tering in a pA collision is larger by a factor A 1//3 than in the corresponding pp collision; 
note that I assume homogeneity in impact parameter space, for simplicity). Thus, the 
ratio © measures the difference between pp collisions and an incoherent superposition 
of proton-nucleon collisions, and is an useful observable to pinpoint collective effects 
like gluon saturation in the wavef unction of the incoming nucleus 3 . 

So, how to compute gluon production in pA collisions ? Leaving aside the fully 
numerical calculations H, [Tjl] (and Refs. therein) based on the classical version of 
the CGC theory (the McLerran-Venugopalan model IlLul]; see also Sect. Ill below), 
which apply to both pA and AA collisions but have not been extended yet to include the 
quantum evolution with r\, all (semi) analytic calculations of the i? p A-ratio 
l22l I3ITI rely on the factorization of the cross-section for particle production in terms of 
the gluon distributions in the target and the projectile ("frr-factorization"), which has 
been since long proven for pp collisions (see, e.g.. lp2TI ^ and has been shown in the 
recent years to also hold for pA collisions lll2Ll33Ll34l 13511 . provided a special definition 
is used for the nuclear gluon distribution [36] which incorporates rescattering effects. 

Up to an irrelevant normalization factor which cancels out in the ratio ©, the gluon 
yield in pA collisions can be computed as (see kinematics in Fig. [3]) 



(5) 



where y denotes the gluon rapidity, related to its longitudinal momentum fraction x via 
y = ln(l/jc), and <PA(k±,y) is the gluon occupation factor 4 , i.e., the number of gluons 
of given spin and color per unit phase-space in a nucleus with atomic number A: 
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Recall that, in pA collisions, we do not expect this ratio to be modified by final-state interactions. 
4 I am glossing here over the subtle (but not essential for the present discussion) distinction between the 
canonical gluon occupation factor (defined as the expectation value of the Fock space number operator 
a\a^ in some specific gauge [T3II '). which is the quantity meant by Eq. (J6J, and the special 'gluon 
distribution' (a scattering amplitude for a color dipole) which strictly speaking enters the fcr-factorized 
expression Q for the gluon production (see, e.g., 1361 13511 for details). 
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FIGURE 6. pA collision at forward rapidity: the (relatively large-*) gluon emitted by the 
incoming proton scatters off the highly-evolved gluon distribution in the nuclear wavefunction 
at small-*. Non-linear effects (which correspond to either gluon recombination in the nuclear 
wavefunction, or multiple scattering of the 'external' gluon) are indicated by small cercles. 



(N c = 3 is the number of 'colors'.) For pp collisions, a similar formula holds with 
<Pa(£l,Ji) replaced by q> p (k ± ,yi). Note that y 1)2 = y ± r], with y = \n(y/s/k±) (cf. 
Fig. |3l). Thus increasing 77 amounts to increasing y\ but decreasing yi. 

Clearly, the nuclear effects enter the previous formulae via the nuclear gluon distri- 
bution <Pa : Unlike its proton counterpart <p p , which is the standard, leading-twist, 'un- 
integrated gluon distribution' generally used in conjunction with ^-factorization IB2I1 . 
q>A encodes higher-twist effect associated with the rescattering suffered by the gluon 
partaking in the collision while it crosses the nucleus. 

At central pseudo-rapidity (ri = 0) and for the high-p^ kinematics at RHIC, y\ = yo is 
rather small, so one can assume that the gluons in the nuclear distribution <PA(k±,y\) are 
produced by direct radiation from the valence quarks. In that case, Eq. © describes the 
multiple scattering between a gluon in the proton and the valence quarks in the nucleus, 
as illustrated in Fig. El The effect of the rescattering on the ratio R p a will be discussed in 
the next section, within the context of a simple model for the distribution of the valence 
quarks, introduced by McLerran and Venugopalan yj]. 

At forward, and sufficiently large, pseudo-rapidities (77 = 2 — 3), the gluons with yi = 
yo + rj in the nuclear wavefunction are predominantly produced by quantum evolution 
from lower values of yi 113 711 . that is, they are emitted by gluons with larger values of x 
which are themselves radiated by the valence quarks. Since the gluon density was rather 
large to start with (due to the large number of valence quarks), and it is further amplified 
by the evolution with Tj, the newly emitted gluons propagate in a strong background 
field — the color field created by sources (valence quarks or gluons) at larger values 
of x — and undergo multiple scattering off the latter (see Fig. |6j). This rescattering 
introduces a negative feedback (it reduces the number of gluons), and eventually leads 
to gluon saturation. Computing in this non-linear regime requires an analysis of 
the quantum evolution in t he p resence of strong fields, which has been completed only 
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with results to be discussed in Sect. IV. 



The simplify the discussion while keeping all the salient physical features, it is 
convenient to replace the ratio in Eqs. ©-© by the simpler quantity 



% pA[ k L ,y) = J,f)f (7) 
A l P<p p (k±,y) 

to be considered for y = y%, i.e., for the rapidity of the gluon in the nuclear wavefunction 
which participates in the collision. Eq. measures the difference between the gluon 
distribution in the nucleus and that in the proton (scaled up by A 1 / 3 ) at the same (small) 
value of x, and thus is the most direct expression of the nuclear effects which should 
be responsible for the deviation of the experimentally measured ratio R p a from one. In 
particular, qualitative effects like the Cronin enhancement or the high-p^ suppression 
should be already visible (and theoretically easier to study) on the ratio ©; this is 
confirmed by the analyses in Refs. lEol l22ll . which have found a similar qualitative 
behaviour for the quantities defined in Eqs. ©-© and, respectively, Eq. ©. In what 
follows I shall exclusively focus on the simpler ratio, Eq. 0. 



III. CENTRAL RAPIDITY: CRONIN PEAK IN THE MV MODEL 

So long as the rapidity y is not too large, one can ignore quantum evolution towards small 
x and describe the gluon distribution in the nucleus as the result of classical radiation 
from the valence quarks. This is probably a good approximation for particle production 
in d+Au collisions at central rapidity (rj = 0), since in that case y = yo ~ 3 is still 
quite small. (From experience with the phenomenology at HERA, we expect small-* 
evolution effects to become important only for y > 5, corresponding to x < 10~ 2 ; see 
e.g. 114211 and references therein.) 

The gluon spectrum dNA/dyd 2 k±d 2 b± is produced by the valence quarks located 
within a tube of transverse area ~ l/k± which traverses the nucleus in longitudinal direc- 
tion at impact parameter b± (so the length of the tube scales like A 1 / 3 ). If the transverse 
section of this tube is much smaller than the area of a single nucleon (i.e., k 2 ^ ^> Aq CD ), 
then the valence quarks which are included inside the tube belong typically to different 
nucleons. It is then reasonable to assume, following McLerran and Venugopalan (jl, that 
these quarks are uncorrelated with each other. But although the color charges of these 
quarks sum up incoherently, this is generally not so for the corresponding gluon distri- 
butions (unlike it would happen for photon distributions in QED), because the classical 
color fields obey a non-linear equation: the Yang-Mills equation. 

It is convenient to consider first the case of a proton, since there the color fields are still 
weak, and one can rely on the linear approximation. Then, the total gluon distribution is 
simply the incoherent sum of the distributions produced by the iV c = 3 valence quarks, 
which by themselves can be evaluated with the well known formula for bremsstrahlung: 



After multiplying this by N c and dividing the result by the number of (gluonic) color 
states N 2 — 1 and by the proton area kR 2 , one finds the occupation factor as (cf. Eq. ©): 

<Pa(^)-|, th = ^- (9) 

Physically, \x p is the color charge squared of the valence quarks per unit transverse area 
(here, in the proton). Eq. © is valid for k± ^> Q p with Q p ~ Aqcd a non-perturbative 
scale related to confinement. 

Moving to a large nucleus, one finds }1a — A 1 ' 3 /^, but the calculation of the corre- 
sponding gluon distribution is more elaborate, as it requires the exact solution to the 
Yang-Mills equation with a strong color source. This has been solved indeed, with the 
following result EH : 



9a(*l) = d 2 r ± e- ik ^ ^ —= (10) 



1 - exp | - \ a s N c r\ii A In ^ } 
na s N c r\ 

where A ~ Aqcd is a non-perturbative scale introduced by hand to cut off an infrared 
divergence which physically should be removed by confinement. This formula features 
the saturation momentum Q S (A), defined by 

Q 2 (A) = a s N c VA ln^^ ~ A 1 / 3 InA 1 / 3 , (11) 

(note that when r = 2/Q s (A) the exponent in Eq. (flOb is equal to one), as the intrinsic 
transverse momentum scale which separates the linear regime from the non-linear one: 
(i) At high momenta k± ^> Q S (A), one can expand the exponential in Eq. (fTOl) . and 
thus obtain (I only show here the first two terms in this expansion) : 



9 a(*±) * $ + a s N c (!jf) 



2 r k 2 

ln-A+27-2 



for ^>2,(A). (12) 



The first term, which gives the dominant behaviour at high k±, is recognized as the 
bremsstrahlung spectrum, corresponding to independent emissions from the valence 
quarks. This is as expected: At high momenta, non-linear effects become negligible 
because the corresponding modes have only little overlap with each other. Thus, for 
sufficiently high momenta, <Pa(&±) — A l / 3 <p p (k_[_), and the ratio & p A(k±) approaches 
one from the above. 

(ii) At low momenta k± <C Q S (A) (with k± ^> A though), the dominant behavior is 
obtained after neglecting the exponential term in Eq. ( fTOl . and reads 

(p A (k ± ) ~ (ln%^ + ^(1)1 ; for k ± <Q s {A). (13) 
oc s N c { k\ ' ) 



This is parametrically large, &{\/ <X S ), but increases only slowly (logarithmically) when 
decreasing k± and/or increasing the atomic number A. This behaviour, which should 



be contrasted with the corresponding power-law behaviour (in both 1 jk\ and A) of 
Eq. (fT2"l) . is the hallmark of gluon saturation (here, within the MV model). 

Because of this change in behaviour, the ratio © becomes very small at low momenta 
k± <$C Q s (A) (where it is proportional to k 2 ^), but it increases rapidly with z = k 2 J Q 2 (A), 
and becomes parametrically large (and, in particular, larger than one) for momenta in 
the vicinity of the nuclear saturation momentum. Indeed, for z ~ 1, <Pa ~ l/oc s N c and 
<Pp ~ Mp/<2?(A) (cf. Eq. ©), and by also using fi A =A l / 3 /i p , one finds 

& P A(k±) ~ Pa > 1 for k ± ~ &(A) , (14) 

where (cf. Eq. (fTTT) ') : 

a s NrfiA A 2 

The estimate (fl4l holds up to corrections of <^(1), which are indeed suppressed in 
the limit of a very large nucleus, for which Pa ^> 1. To make contact with RHIC 
phenomenology, it is useful to notice that, for a gold nucleus at RHIC energies, one 
expects Q 2 (A) ~ 2 GeV 2 |5], which together with A ~ 200 MeV implies pA — In 50 ~ 4. 

The previous considerations show that the ration M pA must have a maximum as a 
function of k±, with the position of the maximum near Q S (A) and the height of the 
peak of ^(lnA 1 / 3 ). This is confirmed by an exact (analytic) calculation Bill of the 
nuclear gluon spectrum in the MV model, Eq. (fTUI) . which together with Eq. © for the 
corresponding spectrum in the proton leads to the ratio M pA displayed in Fig. [7] (the left 
figure there). The location zo of the maximum and its height & max (A) = & p a(zo) can be 
computed analytically, in an expansion in powers of 1 / pA 0"lll : 

882 

zo = 0.435 + + ^(p A 2 ), ^max(A) = 0.281 p A + 0.300 + e{p A l ) . (16) 

Pa 

In Fig. |7J it is also shown the result of a calculation including running coupling ef- 
fects in the MV model 5 (the right figure there). This is obtained by simply replacing the 
explicit factors of (X s in Eq. (fTOb by the one-loop running coupling : oc s — > a s (4/r^), 
with <X S (Q 2 ) = bo/ln(Q 2 /AQ CD ). Then, the scale which plays the role of a saturation 

momentum reads simply Q 2 (A) = boN c HA, and scales like A 1 / 3 . In that case too, the 
gluon distribution in the MV model can be computed analytically and the results 
are rather similar (at least, qualitatively) to the fixed coupling case. In particular, the pa- 
rameter pA = In Q 2 (A)/A 2 plays again an important role — with a running coupling, this 
appears even more naturally, as the inverse of the running coupling: pA = bo/a s (Q 2 (A)) 
— , and expansions like those in Eq. (fToT) can be again derived. (It turns out that the co- 
efficients of the dominant terms in these expansions — for both zq and M max (A) — are 
the same for fixed and running coupling.) 



The inclusion of running coupling effects in an otherwise classical calculation may look as rather ad- 
hoc. Still, this has the merit to approximately include an important class of quantum corrections, which 
are potentially large. 




FIGURE 7. The ratio & p a(z) as a function of the scaled momentum variable z = k 2 / Q 2 (A) in 
the fixed (left) and running (right) coupling McLerran-Venugopalan model for Pa = 6. The thick 
line corresponds the ratio 3% p a{z)\ the solid line shows the saturation contribution M p ^[z)\ the 

dotted line shows the twist contribution ^!™ st (z). (From Ref. [31].) 



The plots in Fig. |7] not only confirm that a maximum exists, but also show that this 
maximum is rather well pronounced, which could not have been anticipated solely on the 
basis of the previous approximations, Eqs. (fT2t and (fT3t . To understand that, one needs 
also the behaviour of the nuclear gluon spectrum <Pa(&±) for momenta around Q S (A). 
The calculations in Refs. J^HS] reveal that, when increasing k± above Q S (A), q>A{k±) 
starts by decreasing exponentially with z = k 2 /Q^A), before eventually relaxing, for 
z > In Pa, to the power law decay displayed in Eq. (IT2h . To better illustrate this, Fig. |7] 
exhibits also the individual contributions to & p a denoted as ffl*^ and ^^ lst which are 

obtained after separating <Pa(&±) into two pieces: the 'twist' piece (p^ lst , which resums 
all the terms which at high-momenta decay as inverse powers of z (the first two terms 
in this series are shown in Eq. (TT21) ). and the 'saturating' piece <p| at = <Pa — <p^ WISt , which 
at z <C 1 has the logarithmic behaviour shown in Eq. (fT3t . whereas for z > 1 decays 
exponentially with z. For any z < lnp^, the saturating piece dominates over the twist 
piece, and thus M p a has a rapid fall off at momenta just above the maximum (within the 
range 1 ^ z ^ In Pa), leading to the well-pronounced Cronin peak manifest in Fig.|7] 

To conclude this discussion, let me comment on the physical origin of the Cronin 
enhancement in this classical model for saturation. As already mentioned in the Intro- 
duction, the Cronin peak is generally associated with multiple scattering within the nu- 
cleus (cf. Fig. El). This is also the case in the present context, except for the fact that, 
since we are looking directly at the wavefunction (rather than at a scattering process), 
the multiple scattering is replaced by non-linear effects in the production of virtual glu- 
ons from the valence quarks. (One can say that, after one gluon is produced by some 
quark, it scatters off the color fields produced by the other quarks.) Since the gluon 
mutual interactions are repulsive, and are stronger at low momenta, their net effect is 
to rearrange the gluon distribution in transverse momentum space 112(1 13 ill: Whereas in 



the bremsstrahlung spectrum oc \jk\ most quanta are located at low momenta, in the 
presence of non-linear effects, the would-be low-^ gluons are pushed towards higher 
momenta, in such a way to minimize their repulsion. Some of these 'displaced' glu- 
ons are responsible for the 'higher-twist' contributions to the tail of the distribution at 
k± 3> Q S (A) (cf. Eq. fiY2\i ). But most of them are quasi-uniformly redistributed at lower 
momenta ^ QsiA), thus giving a saturation plateau at k± < Q S (A) (cf. Eq. fil3\> ) 113 ill . 



IV. NON-LINEAR GLUON EVOLUTION IN THE CGC 

The non-linear evolution of the gluon distribution with increasing y = In I fx is described 
by the renormalization group equation for the CGC, also known as the JIMWLK equa- 
tion which is a functional (or operator) equation, that is, it generates an 
infinite hierarchy of ordinary integro-differential equations which couple the evolution 
of various n-point functions. These equations resum those radiative corrections which 
are enhanced by either the large rapidity gap (namely, the terms of order (oc s y) n for 
any n > 1), or the high gluon density (since, e.g., the gluon occupation factor being of 
&(l/oc s ) at saturation, it interferes with the perturbative expansion). This complicated 
evolution is illustrated in Fig. |5] which shows a typical gluon cascade which develops in 
the nuclear wavefunction at small x. The horizontal rungs in this cascade are radiated 
gluons which are strongly ordered in rapidity, while the mergings between various verti- 
cal branches, as well as the multiple scattering of the produced gluon, are representative 
for the non-linear effects in this high-energy environment. The general solution to the 
JIMWLK equation is not known, but approximate solutions have been constructed which 
separately cover the non-linear regime deeply at saturation, k± <C Q s {A,y), (where re- 
markable simplifications occur in the evolution equations, due to saturation [40]), and 
the linear regime at k± ^> Q s (A,y), where the gluon occupation factor <PA(k±,y) obeys 



a closed, linear, equation: the BFKL equation [37]. 

Remarkably, it turns out that the gluon occupation factor deeply at saturation (k± <C 
Q s (A,y)) has the same simple form as in the classical MV model (cf. Eq. (fT3l ). namely, 

(p A (k±,y) « ^ { ln 6j ^ ,y) + ^C 1 )} » for k ± <^Q s (y,A), (17) 

where the saturation momentum is exponentially increasing with y [Q] : 

Q]{A,y) ~ Q 2 S (A) e^ , (fixed coupling), (18) 

where c ~ 4.88 and d s = a s N c /n. In this and the subsequent formulae, the variable y 
denotes the difference from the rapidity yo corresponding to r\ = 0; that is, from now on, 
y is numerically the same as the pseudo-rapidity r\ . 

Note the special way how Eq. (fTTT) depends upon A, y and k± : <PA(k±,y) is solely a 
function of the dimensionless ratio z = k\l Q%(A,y). This property, known as geometric 
scaling ll25ll . reflects the fact that Q s (A,y) is the only intrinsic scale at saturation. 

Another remarkable feature is that geometric scaling is approximately preserved 
within a wide range of momenta above Q s , where the evolution is linear II 1 3fl . One finds 



indeed that, within the range 

Q s (A,y) < k± < Q 2 (A,y)/Q s (A) ('extended scaling window'), (19) 
the solution to the BFKL equation with saturation boundary conditions at k± ~ Q s is 



well approximated by the scaling form H13LI141I24I1 : 



^ )K i{ ln efc) + ^)}(T) r ' ™- <20) 

The difference 1 — 7 « 0.37 is sometimes referred to as an 'anomalous dimension'. 
Indeed, the fact that 7 is strictly smaller than one makes the function in Eq. (OUb to 
show weaker dependencies upon A and \/k\ than the leading twist approximation (that 
one could naively expect to apply at momenta above Q s ). Below, I shall succinctly refer 
to the window (fT9l) as the 'BFKL regime' . 

The standard leading twist (or DGLAP) approximation is recovered only at even 
higher momenta k± ^> Q~(A,y) / Q S {A), where the gluon occupation factor is extremely 
small, and the only trace of saturation is visible in the fact that it is the original saturation 
momentum at y = 0, i.e., Q S {A), which acts as the infrared cutoff for the transverse 
phase-space available for evolution. Namely, because of the collinear singularity of 
QCD (see, e.g., Eq. ©), the phase-space for the emission of a single (high-fc^) gluon is 
dkj^/k 2 ^, and therefore the total phase-space available for the evolution from Q S (A) up 
to k± is p(A,k±) = ln£j_/<2^(A). Within the present approximations, the solution to the 
linear equation which resums such small-* and collinear gluon emissions is given by 

<PA(k ± ,y) ~ ^exp{y/4a s yp(A,k ± )} , (21) 

which is recognized as the evolution of the leading-twist term in the classical distribu- 
tion (the first term in the r.h.s. of Eq. ([721 ). Eq. (OTI) is known as the "double-logarithmic 
accuracy", or DLA, approximation, and is strictly valid so long as a s yp 3> 1. 

By using the previous formulae, together with similar formulae for the proton, it is 
possible to compute the ratio & P A(k±,y) and follow its evolution with y. The corre- 
sponding formulae for the case of a running coupling a s = Ct^k 2 ^) are quite similar 6 (at 
least, within the kinematical ranges in which Eqs. (fTTt and (|2H apply, i.e., at saturation 
and in the BFKL regime), with one crucial difference though: the functional form of the 
saturation momentum, which for running coupling reads lll3Lll4ill6ll 

Q 2 (A,y) ~ A2 CD exp ^Icby + p 2 ^ , (22) 

where b = 127Y C /(1 IN C — 2Nf) and c ~ 4.88 is the same number as in Eq. (fT8l . As 
compared to the corresponding formula for a fixed coupling, Eq. (fTSl) . the expression 



More complete formulae for the case of a running coupling can be found in the original literature 
H l2l , and are also summarized in Ref . HI . 



(122b shows a less rapid increase with y at high energies, and also a weaker dependence 
upon the atomic number A, which becomes less important with increasing y [16]: 

Q%(A,y) \ P 2 A-P 2 P \ , . , ^ 9 ^ Jl „~ 

WFJ)^ QXV \um) for 2cb y»p*>p P > ^ 

and eventually disappears: 

(g(A,y) ~ A^ CD e^ for 2dy » p* . (24) 

For sufficiently large energies, quantum evolution with running coupling washes out 
completely the difference between a nucleus and a proton ! 



V. FORWARD RAPIDITIES: HIGH-p ± SUPPRESSION 

We are now prepared for a study of the evolution of the ratio & P A(k±,y) with increasing 
y, starting with the initial condition provided by the MV model (which exhibits Cronin 
enhancement at intermediate momenta, as discussed in Sect. III). 



V.l. Quantum evolution of M p a : General features 

Let me first summarize the main features of the evolution, and then develop some of 
them in the next subsections: 

a) The main effect of the evolution is a rapid suppression of the ratio & p a, due to 
the different evolution rates for the gluon distributions in the nucleus (the numerator in 
Eq. 0) and in the proton (the denominator there). The proton distribution grows faster 
because, for the same values of k± and y, the transverse phase-space available for its 
evolution is larger than that for the nucleus. Indeed, as noticed in Sect. IV, the transverse 
phase-space f k± (d 2 p^/p^_) is limited by the infrared cutoff introduced by the initial 
conditions at y = 0. For the nucleus, this cutoff is the relatively hard scale Q S (A) associ- 
ated with saturation ; thus, the nuclear phase-space p(A,k±) = \nk\/Q%(A) is consid- 
erably smaller than the proton one, p(/?, k±) = lnfc 2 /Q 2 with Q p ~ Aqcd- Specifically, 

k 2 k 2 2 (A) 

p(p,k ± ) = ln-£ = In-^ + ln^^ ~ p(A,* ± )+p A> (25) 
sd p Us (/v U p 

where p^ ^> 1 has been generated according to Eq. (fT5t . Correspondingly, M P A(k±,y) 
decreases very fast with y, and already after a short evolution 7 y ~ 1/(6c s Pa) ~ 1 it 
becomes smaller than one at all but the asymptotic momenta. 

b) By the same argument, the suppression goes away at extremely large momenta, 
where the difference between Q S (A) and Q p becomes unimportant in computing the 



7 Recall that a s {Q](A)) = b/p A - 1 /p A since b = &{\). 



phase-space. In fact, when k± ^> Q s (A,y), one can use the DLA fonnula (OTT) for both 
the proton and the nucleus, and thus deduce (for fixed coupling) : 



m pA {k^y) ~ e-'WWW for k ± ^Q s (A,y), (26) 

which approaches one from below when k± — > °°. 

c) The suppression rate dln^^/dy is largest at small y and for not so large trans- 
verse momenta [say, for k± ^ Q%(A,y) / Q S (A)], since in this regime the dissymmetry 
between the evolution of the proton and that of the nucleus is most pronounced: The 
proton is in the DLA regime, and thus evolves very fast (because of the large transverse 
phase-space available to it), whereas the nucleus shows geometric scaling, and evolves 
only slowly (because, so long as a s y < 1, the nuclear saturation momentum rises very 
little; see Eq. (fT8K 

This explains, in particular, the rapid suppression in M p a observed in the early stages 
of the evolution in the numerical study in Ref. Il22ll . 

d) For larger y such that a s y ^ 1, the ratio & P A{k^y) is monotonously increasing 
with k±. That is, the Cronin peak has flattened out during the first 1 / a s units of rapidity. 

e) The flattening of the Cronin peak cannot be attributed to the proton evolution alone 
— the latter produces a quasi-uniform suppression in M p a at momenta around Q s (A,y), 
so, by itself, it would preserve a local structure like a peak — , rather this must be related 
to the evolution of the nucleus. As we shall see, it is indeed the nuclear evolution which 
washes out that distinguished feature of the initial distribution which was responsible for 
the existence of a well-pronounced peak at y = : the exponential fall off of the gluon 
occupation factor at momenta just above the saturation plateau (cf. Sect. III). 

f) Whereas the generic features of the evolution, as described above, are qualitatively 
similar for both fixed and running coupling, important differences persist between these 
two scenarios as far as the details of the evolution, and also the precise structure of the 
final results, are concerned. Specifically, after including running coupling effects, the 
evolution appears to be slower (one needs a larger increase in rapidity to achieve a given 
suppression in M pA ), but eventually stronger (the final value for M p a which is obtained 
after a very large evolution in y is significantly smaller with running coupling than with 
fixed coupling). Let me be more specific on these two points: 

To appreciate how fast is the evolution, let me introduce the rapidity yo after which the 
ratio M p A(k±,y) at k± ~ Q s (A,y) decreases from its initial value of &{pA) (cf. Eq. (O) 
to a value of &{\). In the next subsection, we shall see that 

1 ln 2 p A (lnlnA 1 ^) 2 

y°~JW ~^ i 7T7s — ( fixed coupling) , (27) 

4a s Pa In A 1 / 5 

(which incidentally is a very small rapidity interval: a s yo <C 1), and, respectively, 

yo — yr hi Pa ~ hi A 1 7/3 (running coupling), (28) 
4b 

which for large A is parametrically larger than the fixed coupling estimate in Eq. d27l) : 
thus, the running of the coupling slows down the evolution. 



Furthermore, to characterize the strength of the suppression after a very large rapidity 
evolution, consider the limit of M p a when y — > °° with fixed z = k\ /Q 2 s {A,y) . (This is the 
meaningful way to take the large-j limit, since the interesting physics is located around 
the nuclear saturation momentum.) For z= 6{S), one finds: 

& P A (z ~ 1 , y -> oo) ~ ( A i/3 p A y-y ( fixed cou P lin s) > ( 29 ) 
and, respectively, 

& P A (z ~ 1 , y -> 00 ) = -^3 (running coupling) . (30) 

As anticipated, for large A, the running coupling result (l30l is much smaller than the 
corresponding one for fixed coupling, Eq. (|29b (recall that 1 — 7 ~ 0.37). 

In fact, the power of A 1 / 3 in the r.h.s. of Eq. QUI) is simply the factor introduced 
by hand in the definition of M p a- That is, the result (137)1) arises directly from the 
observation that, with a running coupling and for sufficiently large y, the nuclear and 
proton saturation scales coincide with each other, cf. Eq. (EU, so the corresponding 
occupation factors will coincide as well, in the whole kinematic range for geometric 
scaling (which includes the saturation domain at z < 1 and the BFKL regime (fT9l)). 

g) The dependence of the ratio M p a upon A is also interesting, since this corresponds 
to the centrality dependence of the ratio R^au measured at RHIC lEM Q]. Consider 
the A-dependence for momenta around the Cronin peak: Whereas at y = 0, the ratio 
&pA(k± ~ Qs(A)) is logarithmically increasing with A (recall Eq. (O), this tendency is 
rapidly reversed by the evolution (see the next subsection) : After only a small rapidity 
increase y ~ \/(u s Pa) ~ 1> &pA (^_Ljj) becomes a decreasing function of A for any 
k±, in qualitative agreement with the corresponding change in the centrality dependence 
observed in the data, cf. Fig.|5j 



V.2. The suppression of the Cronin peak 

I shall now use some simple calculations to illustrate the prominent role played by 
the proton evolution for the suppression of M p a (especially at small y). For definiteness, 
I shall consider the fixed coupling case and focus on transverse momenta of the order 
of the nuclear saturation momentum, since this is the region where the Cronin peak is 
located at y = 0. (The generalization of the discussion below to a running coupling and to 
arbitrary transverse momenta can be found in Ref. IBIID So, let's consider the evolution 
with y of the following quantity (the ratio along the nuclear saturation line) : 

® 9 *{A,y)=& pA .{k L = Q s {A,y),y). (31) 

Along this line, the nuclear gluon distribution reduces to a constant, due to geometric 
scaling (cf. Eq. ([171) ) : 



<PA(k ± = Q s (A,y),y) ~ l/a s N c . 



(32) 



As for the corresponding distribution in the proton, this is given by the linear (BFKL) 
evolution, since k± = Q s (A,y) ^> Q s (p,y)- Specifically, for relatively small y, such that 
Otsy <^ Pa, the evolution is dominated by the large transverse logarithm p(p,k±) ~ pA 
(cf. Eq. d25t ). and the proton is described by the DLA formula (ED, while for & s y ^ pA 
the proton enters the BFKL regime, where Eq. d20b applies. The 'critical' rapidity 
O s yc ~ Pa at which the proton changes from DLA to BFKL is the upper limit for the 
geometric scaling window ( fT9t . now applied to the proton. 

I) a s y <C pA '■ The early stages of the evolution (proton at DLA) 

Using Eq. dZTt with A — > p and k± = Q s (A,y), together with p(p, Q s (A,y)) ~ pA (cf. 
Eq. (1251) ) and the definition (fT31) . one immediately finds: 



Fory = 0, this is parametrically large, M sat (A, y = 0) ~ pA, as expected (cf. Eq. (O). But 
when increasing y (with ct s y <C pA though), the ratio decreases very fast — the DGLAP 
increase of the proton distribution being faster than the BFKL increase of the nuclear 
saturation momentum, cf. Eq. (fTSt — , and becomes parametrically of &{\) already 
after the very short rapidity evolution 



which is Eq. (1271 . This value yo is so small that one can in fact ignore the corresponding 
evolution of the nucleus: The rapid decrease in the height of the peak in the very early 
stages of the evolution is entirely due to the DGLAP evolution of the proton. 

In fact, Eq. (EHt shows that, larger is p A (i.e., higher was the original peak at y = 0), 
faster is the suppression seen when increasing y (i.e., smaller is yo). This is so since pA 
also fixes the transverse phase-space for the DGLAP evolution of the proton, and as 
such it enters the exponential factor in Eq. d55t . For the same reason, the ratio (133b turns 
rapidly into a decreasing function of pA (and thus of A), as anticipated in Sect. V.l. 

I) cc s y <i pA : Proton in the BFKL regime 

For y > y c , the proton enters the scaling window ( fT9t . where Eq. (l20b becomes 
appropriate. I shall shortly argue that, when this happens, the Cronin peak has already 
disappeared; but it is still interesting to follow the ratio & P A(k±,y) further up along 
the nuclear saturation line. By using Eq. (l20t with Q s (A,y) — > Q s (p,y), together with 
Eq. (l32t and the relations (which follow from Eqs. (fTTT) . (fT31) and Ha = A 1//3 jU p ) 




(33) 



a s yo 



In 2 Pa 
4pA 



« 1, 



(34) 



a 2 (A,y) _ Q?(A) 




In 



Q 2 s(A) 



Pa, 



(35) 



Q 2 s {p,y) Qj 




one easily finds 



as anticipated in Eq. Note the power of A in the denominator: this provides a strong 
suppression factor which is independent of y. This power was missing in DLA (compare 
to Eq. (1331 ). but appears here as a consequence of the 'anomalous dimension' 1 — 7 > 
characteristic of the BFKL solution in the vicinity of the saturation line lfl3lH4ll . 

Note furthermore that the result in Eq. d36t is independent of y : the y dependencies 
have cancelled in the ratio (l35l) between the nuclear and the proton saturation momenta. 
Thus, as compared to the (proton) DLA regime at y < y c , where & sa t(A,y) is rapidly 
decreasing with y, in the BFKL regime at y > y c this ratio stabilizes at a very small 
value, proportional to an inverse power of A. 

At this level, one can easily anticipate that the behaviour at large y should be quite 
different with a running coupling: Indeed, in that case, the y-dependencies do not 
compensate in the ratio Q 2 s {A,y) / Q 2 s (p,y) (cf. Eqs. dZ2b-dZ5b'). so the function ^ sat (A,y) 
keeps decreasing also in the (proton) BFKL regime, down to the limiting value ^ sa t = 
1 /A 1 / 3 , cf. Eq. (l30l) . In practice, this value is reached for 2cby > p\ (cf. Eq. d23h). 



V.3. The flattening of the Cronin peak 

To study the fate of the Cronin peak with increasing y, one needs to extend the 
previous analysis to momenta outside (but near to) the saturation line k± = Q s (A,y). 
In this region, the rapid evolution of the proton provides a strong suppression, as already 
discussed, but by itself this evolution is quite uniform in k± and it could preserve 
a local peak (whose height would be rapidly decreasing though). From Fig. |U we 
note that the current data are not accurate enough to tell us whether the peak actually 
persists with increasing y, or not. But from a theoretical perspective we expect this 
peak to rapidly flatten out and disappear after only a short evolution, for reasons to 
be explained shortly Such a behaviour has been first seen numerically, in Ref. 
1E2I1 . where the BK equation I29I l30ll has been solved with initial conditions of the 
MV type (cf. Sect. III). This scenario is further supported by the analytic estimates in 
Sect. IV, which show that, when extrapolated down to k± ~ Qs(A,y), the BFKL solution 
(EH becomes parametrically large, of &{\/ cc s ), and thus can be directly matched onto 
the corresponding extrapolation of the solution at saturation, Eq. (ITTb . This property 
suggests that there is no room left for a (parametrically enhanced) peak around Q s (A,y). 

To better appreciate why this matching issue is relevant for the problem of the flat- 
tening of the Cronin peak, one should recall from Sect. Ill that the existence of a well 
pronounced peak at y = is precisely related to the large mismatch, around k± = Q S (A), 
between the 'saturating' distribution <p| at and the 'twist' distribution <p^ wlst (the sum of 
the power-law tails at high momenta) : For k± ~ Q S (A), <p| at ~ l/a s is parametrically 
larger than <p^ wlst ~ l/(a s Pa) ~ 1 (cf. Eqs. (fT2t and (fT3t). Because of this mismatch, 
there is a sudden drop (actually, an exponential fall-off) in the nuclear gluon distribu- 
tion from the saturation plateau at low k± to the power-law tail at high k±. This rapid 
fall-off is responsible for the well pronounced peak in the ratio M v a visible in Fig.|7J 

One can now understand why the quantum evolution leads to the flattening of the 
Cronin peak: Since non-local in k^, the evolution replaces any exponential profile in the 
initial conditions (here, the one at momenta just above Q s ) by a power-law profile. Thus, 
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FIGURE 8. The Cronin ratio M p a{z) (with z = k\ /Q 2 S (A)) below and near the saturation scale 
for Pa = 6. The solid lines correspond to an evolved nuclear wavefunction by AY <C 1 . The dotted 
lines correspond to an unevolved one (MV). The proton wavefunction is always given by the 
DLA solution. The curves, from top to bottom, correspond to AY = n/(2p^) with n = 0,1,. ..,4. 



with increasing y, the 'exponential gap' between the saturation plateau and the power- 
law tail is rapidly filled up, predominantly due to radiation from those gluons which 
were originally at saturation. Correspondingly, the peak flattens out, and completely 
disappears after a (rather short) rapidity evolution a s y ~ 1 113 ill . 

To explicitly follow the flattening of the peak, one cannot rely on the approximate 
solutions presented in Sect. IV — the latter apply only for larger y with a s y > 1, and 
merely show that the Cronin peak has already disappeared by then — , but rather one 
needs more accurate solutions. In the analysis in Ref. 113 ill , the BK equation has been 
iterated exactly (analytically) to follow the first few steps in the evolution of the nucleus. 
By using this evolution together with the DLA approximation for the proton distribution, 
one obtains the results displayed in Fig.[8]for pA = 6 (the length of a step in rapidity is 
taken to be AY = l/(2p^) with Y = 6c s y). For comparison, Fig. [8] also shows the ratio 
which is obtained when the non-evolved, MV model, distribution is used for the nucleus. 
The rapid suppression of the peak, due to the fast rise in the proton distribution, is clearly 
seen in both cases. But in the absence of nuclear evolution the peak is always there; just 
its amplitude gets smaller and smaller. By contrast, when using the properly evolved 
nucleus distribution, the flattening of the peak is manifest, and in fact the maximum has 
almost disappeared already after an evolution AY = 2/ Pa ~ 0.3. 
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